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1. Introduction 

With the on-going heavy ion collision program at RHIC and the forth-coming one at LHC, 
a significant amount of resources is being devoted to the study of hadronic matter at very high 
energy densities. Yet, it is sometimes difficult to express the ultimate physics goals that these 
programs attempt to meet in a conceptually clean manner. After all, nobody expects that one 
would be able to disprove the validity of QCD through these programs, while the computational 
tools that are being developed need to take into account the rather specific background effects and 
non-equilibrium features that may hamper the experimental setup, and are thus somewhat specific 
to these circumstances, without necessarily a much wider applicability in other contexts. 

Of course, it might be that the ultimate motivation for the heavy ion program does not lie 
within QCD phenomenology but elsewhere. In documents directed towards the general public, 
in particular, it is often argued that the main aim is to recreate for a moment the conditions that 
existed in the very Early Universe. In other words, one is hoping to produce a state which is to a 
good approximation thermalised, with a temperature in the range of a few hundred MeV, whereby 
the results should be of interest to the cosmology community. 

It is a fact, though, that most cosmologists have shown relatively little interest in the physics 
related to the QCD epoch since the early 1990's. The main reason is that in order to create cos- 
mological relics, i.e. something that could be astronomically observed today (say, baryonic matter, 
dark matter, inhomogeneities in the various matter distributions), one needs to deviate from thermo- 
dynamical equilibrium. This does not happen easily in cosmology, though, given that the expansion 
rate is very small compared with the rates of microscopic processes. In fact, for strongly interacting 
particles, the only possibility of a deviation from equilibrium is through the existence of a first or- 
der phase transition. But as lattice studies have been indicating since a long time already (for recent 
studies, see refs. [[I]]), the transition in physical QCD seems to be an analytic crossover rather than 
a real phase transition (or, at most, an extremely weak first order transition). Such a scenario does 
not appear to lead to any direct cosmological relics (for a recent review, see ref. [Q]). 

It turns out, however, that this argumentation may be overly pessimistic. In some sense the 
situation is analogous to the case that will be met in the /^-experiments at the LHC: the goal is to 
find "New Physics", be it a Higgs or something more exotic, but QCD still plays an instrumental 
role, because it produces a significant background that needs to be understood extremely well. 

As a first example in this direction, let us mention the properties of the primordial gravitational 
wave background. It is assumed that inflation generates a background with a certain "flat" spec- 
trum. The spectrum observed today depends, however, also on the history of the Universe after 
the inflationary period; in particular it depends on how long a certain wavelength spends inside 
the horizon after re-entering it. This in turn is dictated by the thermal history and the equation-of- 
state of the matter filling the Universe. Consequently, the QCD epoch, during which the expansion 
rate changes by a significant amount, does lead to a distinctive feature in the gravitational wave 
spectrum that could be observed today [||] (if our instruments were precise enough). 

Another example, closer to the topic of this talk, is the problem of Dark Matter (DM). Again, 
the main "signals" are the measurement of the DM relic density and, ultimately, the discovery of 
DM particles. As we will recall in the next section, however, the thermodynamics of the QCD 
epoch does play an important "background" role in determining the relic density. 
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2. The problem of Dark Matter 

A number of independent observations, ranging from the length scales of galaxies (rotation 
curves), through those of galaxy clusters (lensing, large-scale structure), up to cosmological scales 
(anisotropics in the microwave background), have been consistently making a case for the existence 
of non-baryonic DM for quite a while already. The amount of DM is estimated to presently be 
about 20% of the total energy density, with a relative error of at most 5% [|J]. The error gets 
constantly reduced through new observations, and is expected to reach the 1% level through the 
Planck mission at the latest. Of course, all of these observations are indirect, and involve some 
uncontrolled systematic errors. Still, the facts that very different types of observations produce 
consistent numbers, and that the outcome has remained stable for a long time already, make a 
reasonably credible case for the existence of DM. 

It may be amusing to recall the DM relic density in absolute units as well. At the current 
moment in the Universe expansion, it amounts to 

GeV 

m 3 

This may appear to be a small number, but is in fact large compared with the current average energy 
density of the total of all known forms of matter, namely 

GeV 

^Baryon ~ 0.2 5- . (2.2) 

m 3 

Given that all existing evidence for DM is based on its gravitational interactions, the more 
precise nature of DM remains unknown. 1 The most popular candidate is so-called Cold Dark 
Matter (CDM), consisting of particles with a mass m > 10 GeV, related perhaps to supersymmetry 
or supergravity. These particles can be called Weakly Interacting Massive Particles, WIMPs. The 
hope would then be to discover CDM not only through astronomical observation, but also through 
directly creating WIMPs in future collider experiments like the LHC. 

Unfortunately, as of today, there is no concrete evidence in favour of CDM. Therefore the field 
remains open for other candidates as well. In this talk I will concentrate on right-handed "sterile" 
neutrinos as candidates for DM. Within the see-saw scenario it is commonplace to introduce mas- 
sive right-handed neutrinos in any case, to explain the experimental fact that at least two of the 
active neutrinos have a mass. It is then only natural to assume that there are in fact three families of 
right-handed neutrinos in total, like there are of the other particles. It is the lightest among these, 
with a mass in the keV range, which could possible act as DM [jfj[ Because of the small mass, 
and subsequently large average energy at any given temperature, such dark matter is referred to as 
Warm Dark Matter (WDM) in contrast to CDM. 

Now, let us recall briefly why hot QCD does play a role in the determination of the DM 
relic density, even though the DM particles themselves are very weakly interacting (otherwise, 
they would not be "dark"). For CDM, for instance, the relic density is determined through the 
moment when the WIMPs decouple, i.e., when their interaction rate becomes smaller than the rate 

1 In particular, it must be kept in mind the Nature might have chosen an even more exotic explanation for what we 
conceive as DM than just some novel particles: maybe the problems are on the side of gravity, and it is a modification 
of the theory of gravitational interactions which explains the features that we associate with DM. 
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of Universe expansion. The expansion rate, in turn, does depend on the properties of all the particles 
in the plasma, and most of them do feel strong interactions. It turns out, in particular, that WIMPs 
of mass m decouple at T dec ~ m/25 (see, e.g., ref. |01). For m = 10.. .1000 GeV, 7d ec = 0.4... 40 GeV, 
which indeed is a range where quarks and gluons dominate the equation-of-state (EOS). Therefore 
the QCD EOS does affect the CDM relic density on a level which is on par with the observational 
accuracy of the future experiments [§, §]. 

For right-handed WDM neutrinos, the production mechanism is somewhat different from that 
for CDM: these particles do not simply decouple (because their couplings are so weak that they 
never reached equilibrium), but their production rate really peaks at a certain temperature. The tem- 
perature in question depends on the mass of the WDM neutrino, but for the keV range mentioned 
above (see Fig. ^ for the observational constraints which make this range to be the relevant one), 
the peak temperature is around T ~ 200 MeV [||] (see also Sec. |5|). This is just the temperature 
scale of the QCD phase transition or crossover, so that indirect QCD effects on the WDM relic 
density are even more dramatic than on the CDM one. It is these effects that will concern us in the 
following. 



3. Minimal model for right-handed neutrinos 

It has perhaps already become apparent that the philosophies behind the CDM and WDM 
scenarios are rather different. In the CDM case, the hope is that WIMPs exist and offer a window 
to genuinely new physics — the more exotic, the better! In the case of right-handed neutrinos, 
in contrast, the philosophy is to be as down-to-earth as possible: the Minimal Standard Model 
(MSM) is completed by adding only those degrees of freedom which are necessary in any case 
for explaining the established experimental facts concerning neutrino oscillations. Otherwise the 
guiding principle, the construction of the most general renormalizable Lagrangian invariant under 
the gauge symmetry SU(3)cxSU(2)£xU(1)f, remains the same as in the MSM. 

Given these principles, the new Lagrangian reads 

& = ^msm + -M S ]N S - [h as Z a $a R N s + R.c] , (3.1) 

where the generation indices s, a = 1,2,3 are summed over; N s are Majorana fields; M s are real 
Majorana masses; h as are complex Yukawa couplings; L a are lepton doublets; and <j> = is 
the conjugate Higgs doublet. Because of electroweak symmetry breaking, active neutrinos have 
masses in this model. For illustration we will work in the corner of the parameter space where the 
masses are given by the see-saw formula as usual, m Va ~ \h as \ 2 v 2 /M s , where v ~ 246/v2 GeV 
« 174 GeV is the Higgs field vacuum expectation value. 



Now, there is nothing particularly new about the Lagrangian in Eq. (p.l|): it is the very theory 
that is practically always taken as the starting point for the description of massive neutrinos. Usu- 
ally, though, the assumption is made that the Yukawa couplings h as are of the same order of magni- 
tude as the known ones, \h as \ < 1. To reproduce neutrino masses in the range suggested by the os- 
cillation experiments via see-saw, leads then to the assumption that M\, M2, M3 ~ 10 10 ...10 15 GeV. 
Thereby the active neutrino masses are viewed as a window to GUT-scale physics. 

Esthetically, though, a theory with such parameter values may look a bit strange: the (renor- 
malized) Higgs mass parameter is in any case expected to be of the order of the electroweak scale, 
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so one may ask why the Majorana masses should behave any differently. Clearly, it is possible to 
keep the active neutrino masses fixed, if a decrease in the Majorana masses is accompanied by a 
decrease of the neutrino Yukawa couplings in the proportion dictated by the see-saw formula. In 
refs. [0, 0, it was indeed proposed to consider Majorana masses in the range M\ ~ keV, Mi, 



M3 ~ GeV, and Yukawa couplings in the ranges \h a \ \ < 10~ u , \h a 2\, | /i«3 1 < 10~ 7 . Moreover, it 
was pointed out that this possibility leads to some phenomenological benefits compared with the 
usual choice, particularly that the right-handed neutrino Ni now has a lifetime long enough to serve 
as a candidate for DM. The right-handed neutrinos N2, N3 are the ones that induce the observed 
active neutrino mass differences through the see-saw formula, and they also participate in baryo- 
genesis; however, they decay too fast to serve as DM (in fact, they are constrained to decay before 
the Big Bang Nucleosynthesis epoch [|T2|]). Following refs. [10, TTj], we will refer to the extension 
of the MSM with such parameter choices as the "vMSM". 

Given that the Yukawa couplings are small and that consequently the right-handed neutrinos 
interact extremely weakly, we will refer to them as "sterile" neutrinos in the following. For future 
reference, we also define the mixing angles 6 a \ as 



The mixing angles to be considered in the following are all small, \6 a \ | < 10 



-3 



4. Primordial production of right-handed neutrinos 

Even though sterile neutrinos interact very weakly, they do get produced in reactions taking 
place in the Early Universe. This leads to a non-zero primordial abundance, as well as a number of 
potential consequences: 

1. The sterile neutrinos are not stable but decay, for instance through the channel N — ► vy. If 
their life-time is long enough, this results in an X-ray signal which could be observed today. 

2. The sterile neutrinos are massive and thus carry a certain energy density. Provided, again, 
that their lifetime is long enough, this could lead to a contribution to the WDM that plays a 
role in structure formation. 

3. Finally, the various CP-violating scatterings of the sterile neutrinos might contribute towards 
the baryon asymmetry that exists today. 

It is worth stressing that the conventional case of heavy Majorana masses, M\,M2,Mt, ~ 
10 10 ...10 15 GeV, also leads to a certain primordial abundance and, as is well-known, to the possi- 
bility of generating a baryon asymmetry through the mechanism of thermal leptogenesis Jl3|]. This 
region of the parameter space does not lead to the first two consequences, however, since the heavy 
neutrinos decay too fast to still be around today, as sources of X-rays and DM. 

Because of these potential consequences, it is important to determine the abundance of the 
sterile neutrinos as a function of their mass M\ and mixing angles 6 a i, and we now turn to this 
computation. Afterwards, the result of this theoretical computation can be confronted with the 
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observational constraints that follow from the absence of any visible X-ray signal and from known 
properties of large-scale structure (cf. Sec. ||). 

The basic mechanism for the creation of WDM neutrinos is through inelastic scatterings of 
light particles existing in the thermal plasma, followed by an active-sterile neutrino transition. 
These processes were first analysed in detail in ref. [g], and a number of refinements can be found in 
refs. [ 14 ]— [ 18]. Our analysis is based on refs. Jl9| , ^OQ , where the equations entering the production 
rate were derived directly from quantum field theory, allowing for a systematic investigation of all 
the effects that play a role, to all orders in the strong coupling constant. 

It turns out that, as already mentioned and as demonstrated later on explicitly, the production 
rate peaks at temperatures of the order of T ~ 200 MeV <C mw- Therefore, the basic physics that 
plays a role can be understood within the Fermi model. We start by sketching the ideas with a few 
graphs, then discuss the general form of the equations that appear; however, the precise details are 




left out and can be consulted in refs. Q19| , |20| ]. 

An example of a scattering leading to the production of a sterile neutrino N\ is shown by the 
diagram 



(4.1) 



After their production Ni are essentially inert: their average lifetime is longer than the age of 
the Universe, and their density is much below the equilibrium value. Therefore, the N\ produced 
essentially escape the thermal system, just like photons or dilepton pairs produced in heavy ion 
collision experiments do. 

Now, the rate for the production is proportional to the absolute value squared of the amplitude. 
For the process above, this corresponds to the imaginary part of the 2-point function of active 
neutrinos: 



(4.2) 



where the dashed line indicates that the cut, or imaginary part, is to be taken. 

At temperatures of the order of 200 MeV, however, quarks interact strongly, and a perturbative 




evaluation of the 2-loop diagram in Eq. ( |4.2| ) is hardly a good approximation. Fortunately, it is 
possible to express its contents in a more general way, whereby the quark lines can be combined to 



a propagator corresponding to flavour singlet or non-singlet vector or axial current correlators [19]: 




(4.3) 



The vector and axial current parts of this graph can be evaluated also beyond QCD perturbation 
theory, for instance by using chiral perturbation theory or lattice techniques (although the latter are 
faced with the usual problems related to analytic continuation). 
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Finally, there would obviously also be a simpler 1-loop graph that is in principle relevant, 
namely 

w- 

(4.4) 

The imaginary part of this graph is, however, exponentially suppressed by ~ exp(— mw/T), and 
thus insignificant at the temperatures that we are interested in. On the other hand this graph does 
produce a real part (for the real part intermediate particles do not need to be on-shell so that no 
exp(— m\v/T) appears), which indeed plays an important role, as we will see. 




Let us now express these graphs as formulae Q19j |20fl . Denoting by n\ the number density of 
the lightest sterile neutrinos N\, by s(T) the total entropy density, by e(T) the total energy density, 
and by c 2 s {T) the speed of sound squared, the production equation can be written as 



d 

-T — 
AT 



ni(T) 
s(T) 



1 / 3m 2 pi 



3c*(T)s(T) y %ne{T) 



Jd 3 qR(T,q), (4.5) 



where mpi is the Planck mass. The rate, R(T, q), contains in turn the mixing angles and the proper- 
ties of the self-energy L a of active neutrinos of flavour a: 

Here we made use of the fact that sterile neutrinos are on- shell, (q ) 2 - q 2 = M\. Eqs. Q, Q 
show that we need to determine Re£ a , ImZ a , as well as the thermodynamic quantities e(T), s(T), 
and c 2 (r), in order to estimate n\(T). 

Let us start by considering the real and imaginary parts of the active neutrino self-energy. 
The real part originates from the 1-loop graph shown in Eq. ( |4.4J ); a straightforward computation 
followed by an expansion in 1 / leads to the result [ |2l| , [22TJ 

ReE„ = Qa a (Q)+ub a (Q), u=(l,0) (4.7) 



\6GiT 



2 T 4 

F ^0 



b a (Q) = "—q U 20 +cos z w 



T J T V T 



(4.8) 



where is a simple dimensionless function, which can easily be evaluated numerically. The func- 
tion a a (Q) can in fact be ignored, since this part of the self-energy is subdominant compared with 
the tree-level term Q. 

Two interesting remarks can be made related to Eq. (fO|). First of all the parametrically lead- 
ing term, ~ l/ra 2 ^, vanishes. Therefore the result has a prefactor G 2 F ~ gt-/ m w (°f course on ly 
two weak gauge couplings g w appear, whereby we have to divide by after this normalization). 
Second, it turns out that if there are non-zero leptonic chemical potentials in the system, ^ 0, 
then the leading term, ~ l/m 2 ^, no longer vanishes, but produces a term b a {Q) ~ —Gpni, where 



riL is the lepton density [21]. This leads to the possibility of a "pole", or resonance, in Eq. (4.6), 
whereby the production rate can be enhanced by a significant amount [p3|]. Under normal circum- 
stances, however, it is to be expected that the lepton density is of the same order of magnitude as 



the baryon density, ni ~ 10 10 r 3 , in which case these effects are insignificant and can be ignored. 
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As far as the imaginary part of the active neutrino self-energy in Eq. (4.3) is concerned, its 
general structure is 

Im^CGj-G&y^jr^ (4.9) 
Here r is the spatial momentum that flows through the vector or axial current propagator, and J(f 



is a known thermal "kernel", consisting of hyperbolic functions (cf. ref. [19]). The functions p„^ 



represent flavour-singlet and non-singlet vector and axial current spectral functions (i.e., imaginary 
parts of the retarded two-point functions, p = Imn«). 

Now, the spectral functions contain both lep tonic and hadronic contributions. In principle 
these can be computed perturbatively; this leads to structures familiar from the source terms of 



Boltzmann equations. For instance, for the case depicted in Eq. (4.2), 

- «w*huv«i JjS^k /p% * 

x | (2n) i 5 w (P 1 + P2 + P3-Q) n F1 n F2 n F3 A{-m e , m d , -m u ) + Q 

+ {2n) i 8^{P 2 +P3-P1-Q) n F2 n F3 (l - n F1 ) A(m e , m d , -m u ) + 

+ (2tt) 4 5 (4) (Pi + P3-P2-Q) n F m F3 (l - n F2 ) A(-m e , -m d , -m u ) + 

+ {2Tr) 4 S (4 \P 1 + P2-P 3 -Q)n F1 n F 2{l-n F3 )A{-rn, e ,m d ,rn u )+ V 

+ (27r)W\P 1 -P2-P3-Q) ofi(1 - nra)(l - n F3 )A(-m e , -m d ,m u ) + 1— < 



+ (2n) A 5 {i \P 2 -Px-Ps-Q) n F2 (l - npi)(l - n F3 )A(m e ,m d ,m u ) + 
+ (2tv) 4 S^(P 3 -P1-P2-Q) n F3 (l - «fi)(1 - n F2 )A(m e ,-m d , -m u ) + 
+ (27r) 4 5< 4 '(-P 1 - P 2 - P 3 - Q) (1 - n F i)(l - n F2 )(l - n F3 ) .4(fn«, -m d , m u ) J, ^ J 
where 

•A(m e , m<j, m„) = 7 ,J (Pi + m e )7" Tr [(f 2 + m^Y^az^f 3 + m,,)7„aL] . 

Here the labels mean 1 = e, 2 = d, 3 = «. The hadronic contributions should, however, really be 
evaluated beyond perturbation theory, as we have already stressed. 

5. Why does the production rate peak at T ~ 200 MeV? 



Let us recall at this point why the production rate of Eq. ( |4.6[ ) peaks at temperatures of the 
order of a few hundred MeV [0]. 

We have already seen that the real part of the active neutrino self-energy is of the order ReZ a ~ 
GpT 4 a~ l |q| (note that q° ~ |q| for |q| S> M\ as will be the case here). The imaginary part is more 
complicated, but ignoring all masses compared with the temperature, it has dimensionally the form 
Tr[01m^ a ] ~ «p 1 (q°)GjpT 6 f(\q\/T), where / is anon-trivial dimensionless function, numerically 
of order unity. It turns out that the rate peaks at momenta of order |q| ~T, whereby we can replace 
/ by a number of order unity. Thereby the rate becomes 

*(r,q)« y le^i 2 , , 1 F \ . (5.i) 
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It is immediately seen that the result is strongly peaked around temperatures where the two terms 
in the denominator are of similar orders of magnitude, i.e. 

r~(^y~200Me V (-^y . (5.2) 

Thus, for sterile neutrino masses in the keV range, the production rate accidentally coincides with 
more or less the QCD scale. 



6. The role of the QCD equation-of-state 

Having analysed Re£ a and lmL a , we still need to discuss the thermodynamic functions c 2 s {T), 



s(T) and e(T), in order to be able to integrate Eq. (4.5). 

We parameterise the energy density e(T) and the entropy density s(T) through effective num- 
bers of bosonic degrees of freedom, 

K 2 T A 2it 2 T 3 
e(T)= geB (T)— , ^E^r)— . (6.1) 

The energy and entropy densities follow both from the thermodynamic pressure p(T) through 
standard relations, e(T) = Tp'(T) — p{T) and s(T) = p'(T). Furthermore the speed of sound 
squared can be written as 

= ~e r {T) = Tp^Jf) ' (6 " 2) 
In the non-interacting limit, c 2 (T) equals 1/3. To summarise, we need reliable estimates of p(T), 



p'(T) and p"{T) in order to determine the sterile neutrino abundance from Eq. (1.5). 

A reliable determination of p(T) is a long-standing challenge for finite-temperature field the- 
ory. Again, leptonic contributions can be well treated in perturbation theory, while hadronic con- 
tributions, which dominate the structure in p(T) in the temperature range of interest, are in general 
hard to compute precisely. 

Certain limiting values of p(T) are understood better, though. At low temperatures, 
T < 100 MeV, confinement and chiral symmetry breaking guarantee that the system is composed 
of weakly interacting massive hadrons. Treating them as a "gas" of resonances, one can at very 
low temperatures approximate 

rfn-p 4 (£)'•-»■ 

and at somewhat higher temperatures replace this with the corresponding relativistic formulae for 
bosons and fermions, respectively. Even though this prescription is rather phenomenological, lat- 



tice simulations suggest that it may work surprisingly well even up to T ~ 200 MeV [24]. 

For high temperatures on the other hand, T > 1000 MeV, asymptotic freedom guarantees that 
the system can be viewed as a collection of weakly interacting quarks and gluons: 



2^4 



P(T) 



% T 
"90" 



2(Ar c 2 -i)+^\yVf 



[l + ... + 0(g 6 )) , (6.4) 
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Figure 1: Left: g e ff i ^eff as defined in Eq. ( p.l| ), for Nf = 4 QCD with physical quark masses |]27[|. Right: the 
speed of sound squared c%, for the same system. The shaded region is the range of temperatures where our 
recipe is purely phenomenological and needs to be improved through future lattice simulations. 



where N c = 3 is the number of colours and Nf is the number of massless flavours that play a role 
at the temperatures we are interested in. Perturbative corrections are known up to 4-loop order, 
0(g 6 ) [@, [2|], apart from a single missing coefficient. Moreover, quark mass effects, important 
for phenomenological applications, can be incorporated [E7|]. 

In between these two limits the situation is much more complicated. Ideally, one would like 



to make use of lattice simulations of the type in Refs. [|28|]-[]31[]. Unfortunately, it appears that the 
present results are reliable in a fairly narrow temperature range only; for instance, the characteristic 
peak that can be seen in c 2 at around T ~ 70 MeV, due to light pions (cf. Fig. [j](right)), is not visible 
in the existing simulations. In fact the simulations display rather a somewhat deeper dip (down to 
~ 0.1) around the critical region, and then rise at most slightly as the temperature is lowered. 

For these reasons, the results that will be presented in the following adhere to the procedure 



introduced in ref. [27], rather than to lattice simulations. This procedure makes use of a gas of 
hadronic resonances at low temperatures; the most advanced (up to resummed 4-loop level [|25|]) 
weak-coupling results at high temperatures; and an interpolation thereof at intermediate temper- 
atures. 2 Remarkably, the temperature interval where an interpolating function is needed in or- 
der to sew together the two asymptotic functions is fairly narrow, not more than (10. ..20) MeV, 
and centered around T 200 MeV, in curiously good agreement with the crossover temperature 
T c ~ 192 ± 8 MeV as suggested by recent large scale lattice simulations [[32|]. 

The results that follow from this recipe for the quantities that are important for us are illustrated 
in Fig. H for the QCD part, and in Fig. || for the whole Standard Model. In the latter case we have 
displayed, for completeness, the results in a very broad temperature range. 



7. Numerical results 



With the ingredients discussed, Eq. ( |4.5| ) can be solved numerically, as a function of the tem- 
perature. From the result, n\ (T)/s(T), we can derive the current energy density, M\n\ (T), relative 
to the current entropy density. Finally, the current entropy density, whose value is well-known P], 



2 We have corrected a minor error in the numerical results of ref. [f27|] 
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Figure 2: Left: g e ff Aff as defined in Eq. for the Standard Model with m H = 150 GeV Q. Right: the 
speed of sound squared c^, for the same system. The shaded regions are the ranges of temperatures where 
our recipe is purely phenomenological and needs to be improved through future lattice simulations. 



can be traded for the current critical energy density. Thereby we obtain the parameter , char- 
acterising the fraction of the current energy density that is carried by sterile neutrinos. This result 
can then be compared with the observed value for £2 dm- 

In Fig. ^ we show an example of a solution, for specific parameter values [59). We have here 
considered the contribution from the active flavour a = e only. It can be observed that the sterile 
neutrinos indeed get generated at temperatures of a few hundred MeV, and that their relic density 
can be of an order of magnitude which is relevant for the explanation of DM. 

On the other hand, as can be seen in Fig. [3], the current results contain relatively significant 
hadronic uncertainties. These originate primarily from two sources. First of all, as discussed in 
Section |], the imaginary part of the active neutrino self-energy contains the spectral functions of 
vector and axial currents, which involve poorly understood hadronic effects. It is of course clear 
that at low enough temperatures there are no hadronic scatterings taking place, since all hadrons 
are massive. At the same time, at high enough temperatures, the hadronic contributions can be 
computed perturbatively. A very conservative way to estimate the "error" is to consider these as 
limiting cases: as a lower bound no hadrons (N c = 0), as an upper bound (almost) free quarks. In 
the future this error can hopefully be reduced with the help of lattice simulations; the current status 
on estimating some of the relevant spectral functions from the lattice is reviewed in ref. [33]. 

Second, as discussed in Section ^, the hadronic equation-of-state contains significant uncer- 
tainties in the temperature range of interest. It turns out that the most important effect as far as the 
current computation is concerned is the location of the (pseudo)critical temperature T c : whether 
the production rate peaks above or below this temperature has a large effect on the final result, 

I — 3 1/2 

given that the kinetic equation of Eq. Q4.5] ) is inversely proportional to Ih^g^ . We have estimated 
these uncertainties by rescaling the temperature units by 20% in either direction, which certainly 
is a conservative estimate. The current status of lattice determinations of various thermodynamic 
quantities in the vicinity of T c is summarised in ref. [34]. 

The band indicated in Fig. || incorporates both of the error sources discussed, and provides for 
a conservative estimate of the possible hadronic uncertainties in the results. 
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8. Comparison with observational constraints 

We can now confront the theoretical result of Fig. [3] with observational constraints. There are 
observational constraints from two sides. First of all, if sterile neutrinos with a very small mass 
constitute all of DM, then structures on small scales tend to be wiped out, compared with structure 
formation simulations carried out with CDM. Comparing the outcome with actual data, particularly 
in the form of so-called Lyman-a forest observations, which are sensitive to the smallest distance 



scales, puts thus a lower bound on the mass of the sterile neutrinos [p5[]-[|38|]. It appears that the 



lower bound could be as high as (8... 12) keV [2C]. 



Second, the heavier the sterile neutrinos are, and the bigger their mixing angles \6 a \\ with 
active neutrinos, the more likely are they to decay. This leads to a characteristic X-ray signal, which 
however has not been observed. Therefore, it is possible to set an upper bound on a combination 
of the mass and the mixing angles of the sterile neutrinos [JT4 , 39, ff()| ]. 



These two sets of observational constraints are illustrated in Fig. f|, together with the result of 
the theoretical computation outlined above. The basic feature to be observed is that the theoretical 
computation appears to require a larger mixing angle than experimentally allowed, in order to 
produce the observed amount of DM (for «l/T 3 ~ 10~ 10 ). It should be mentioned, however, that 
the structure formation constraints include a number of systematic uncertainties, both on the side 
of the observational Lyman-a data, and for the theoretical reason that they assumed a thermal 
shape for the momentum distribution function of sterile neutrinos, even though these are out-of- 
equilibrium. If the bound happens to be correct despite these uncertainties, Fig. |] indicates that 
sterile neutrinos can only act as dark matter if there is a lepton asymmetry ni/T 3 3> 10~ 10 in 
the Universe, or if there are additional production mechanisms apart from thermal scatterings; it 
could be, for instance, that a certain number density was produced during the inflationary period 
already plj], to which thermal scatterings would then add their contribution. In any case, the band 
around the theoretical curve in Fig. |] indicates the important role that hadronic effects can play. 
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excluded by X-ray constraints 




Mj /keV 



Figure 4: An exclusion plot for the scenario where sterile neutrinos constitute all of dark matter (adapted 
from ref. [SJ; see there also for more details). We compare observational constraints with the theoretical 
computation, with the size of hadronic uncertainties indicated with the band. Here 9 = </£«= e u )T \ 9ai \ 2 - 



9. Conclusions 



The basic scenario of right-handed neutrinos serving as warm dark matter is an old one by 
now It has experienced quite a revival recently, though, because of significant progress both 

on the observational and on the theoretical sides. 

On the observational side, the parameter space (mass, mixing angles) for the lightest right- 
handed neutrino has been strongly constrained by structure formation and X-ray bounds during the 
last year or so. This means that attention can now be focussed on a rather specific situation. 

On the theoretical side, the scenario of right-handed neutrinos serving as dark matter has taken 
a more prominent role, thanks to the realization that there is a minimal model, the "vMSM", which 
not only addresses the dark matter problem, but may also explain neutrino masses [|Kj], baryoge- 
nesis [42j 11], and perhaps also various astrophysical problems p3|]. Furthermore, the theoretical 
tools that are needed in the dark matter computation have reached a more mature level [19]. 

These two sides imply that theoretical computations need to be promoted to a higher level 
of accuracy than before. As has been underlined in this talk, a fully satisfactory analysis in this 
respect is only possible once lattice studies of the QCD equation-of-state and of vector and axial 
current spectral functions with various flavour structures, produce results with controlled statistical 
and systematic errors, such that the band in Fig. |] can be made narrower. Currently the hadronic 
uncertainties are of order 50% (cf. Fig. ||), and it would certainly be desirable to at least half this 
uncertainty. 
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